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ABSTRACT 

Archaeal chromatin proteins share molecular and 
functional similarities with both bacterial and eu- 
karyotic chromatin proteins. These proteins play 
an important role in functionally organizing the 
genomic DNA into a compact nucleoid. Cren7 and 
Sul7 are two crenarchaeal nucleoid-associated 
proteins, which are structurally homologous, but 
not conserved at the sequence level. Co-crystal 
structures have shown that these two proteins 
induce a sharp bend on binding to DNA. In this 
study, we have investigated the architectural 
properties of these proteins using atomic force mi- 
croscopy, molecular dynamics simulations and 
magnetic tweezers. We demonstrate that Cren7 
and Sul7 both compact DNA molecules to a similar 
extent. Using a theoretical model, we quantify the 
number of individual proteins bound to the DNA as 
a function of protein concentration and show that 
forces up to 3.5 pN do not affect this binding. 
Moreover, we investigate the flexibility of the 
bending angle induced by Cren7 and Sul7 and 
show that the protein-DNA complexes differ in flexi- 
bility from analogous bacterial and eukaryotic 
DNA-bending proteins. 

INTRODUCTION 

Organisms in all three domains of life need to compact 
and functionally organize their genomic DNA into the 
relatively small volume of a nucleus or a cell. 
Architectural proteins (histones and other chromatin 
proteins in eukaryotes and nucleoid-associated proteins 



in bacteria and archaea) play an important role both in 
compaction and functional organization of DNA, thus 
affecting DNA transactions as diverse as transcription, 
repair and replication (1). 

In eukaryotes, DNA is wrapped around histone 
octamers, forming nucleosomes. With the aid of other 
chromatin proteins, fibres with nucleosomes are folded 
into higher-order structures, obtaining multiple levels of 
organization. Bacteria organize their genomic DNA into a 
nucleoid, which is shaped by the action of small chromatin 
proteins that bend or bridge the DNA [e.g. histone like 
protein from E. coli strain U93 (HU), integration host 
factor (IHF), factor for inversion stimulation (FIS) and 
histone-like nucleoid structuring protein (H-NS)] (2,3). 
Archaea, constituting the third domain of life, synthesize 
numerous nucleoid-associated proteins, with molecular 
and functional similarities to both bacterial and eukary- 
otic chromatin proteins (4,5). The two main archaeal 
phyla, Euryarchaea and Crenarchaea, express different 
sets of chromatin proteins, none of which is conserved 
throughout the whole archaeal domain. Euryarchaea syn- 
thesize true tetrameric histone homologues, which form 
nucleosomes similar to eukaryotic tetrasomes (6,7). 
Crenarchaea, on the other hand, do not synthesize 
histone homologues, although there are some rare excep- 
tions (8). Each crenarchaeal species encodes at least two 
different small chromatin proteins and several paralogues, 
which may act in concert to compact the genomic DNA 
and to regulate its accessibility (9). 

In the crenarchaeum Sulfolobus solfataricus, four 
nucleoid-associated proteins have been identified so far, 
Alba, SsolOa, Cren7 and Sso7d (the Sul7 family of chro- 
matin proteins). Alba forms dimers in solution (10) and 
has been shown to bridge DNA at low concentrations, 
providing a means of organizing and compacting DNA 
(11) (Laurens el al., submitted). SsolOa homologues 
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(12,13) also exist as dimers in solution and have been 
shown to form protein-DNA complexes, similar in 
appearance to Alba-DNA complexes (11). Otherwise, 
little is known about the role of SsolOa in compacting 
and organizing genomic DNA. 

Cren7 and Sul7 are both small (~7kDa) basic mono- 
meric proteins, which bind to DNA with no apparent 
sequence specificity. Although they share no similarities 
at the amino acid level, their tertiary structures and their 
biochemical properties are similar (14-20). Both proteins 
are folded similarly, containing two antiparallel (3-sheets. 
They mainly differ in the presence of an extended loop 
located in between the two P-sheets in Cren7 (19) and an 
additional C-terminal a-helix in Sul7 (21). Protein-DNA 
co-crystal structures showed that they bind to DNA by 
intercalation of hydrophobic side chains into the minor 
groove (see inset Figure 1), which results in bending of 
DNA by up to 50-60° (23-27). 

The transcriptome of 5. solfataricus P2 shows high 
levels of transcription of both Cren7 and Sul7 (28), and, 
consequently, both proteins are abundant in the cell [~ 1 % 
and up to 5% of the cellular protein for Cren7 (19) and 
Sul7 (20) respectively]. Native Cren7 and Sul7 proteins are 
post-translationally modified by methylation of several 
lysine residues. Nevertheless, the function of methylation 
remains unclear, as it does not change the binding affinity 
to DNA in vitro (15,19), in contrast with the acetylation of 
Alba (29). 

Proteins that bend DNA by sequence non-specific 
binding into the minor groove play an important role in 
DNA compaction and chromatin organization through- 
out all domains of life (1,30-33). In eukaryotes, high- 
mobility group (HMG) proteins form moderately flexible 
bends in DNA by intercalating into the minor groove (31), 
enhancing the overall flexibility of the DNA, and they are 
likely to play a supplementary role in chromatin organiza- 
tion (30). In bacteria, DNA-bending proteins are crucial 
players in chromatin organization, as bacteria lack histone 
homologues. 

Whether Cren7 and Sul7 function similarly to this 
group of bacterial and eukaryotic DNA-bending 
proteins, as suggested by the co-crystal structures, has 
not yet been investigated. A functional difference 
between the two proteins is expected, as they coexist in 



Cren7 Sul7 




A/. 



. £ 

Figure 1. Schematic illustrating the reduction of end-to-end distance 
because of binding of a DNA-bending protein. Binding of a protein 
that induces a bend of deflection angle a reduces the end-to-end 
distance with Ar = — ■ C(ot). Inset: average conformation of the 
Cren7-DNA and Sul7-DNA complexes obtained from MD simula- 
tions. Both proteins bend the DNA by intercalating in the minor 
groove. Images are generated using visual molecular dynamics 
(VMD) software (22). 



species from the order Sulfolobales, but previously 
reported biochemical assays have not yet revealed a sig- 
nificant difference in protein-DNA interactions between 
Cren7 and Sul7. To obtain a detailed understanding of the 
architectural properties of the two proteins, we use a set of 
different techniques, including atomic force microscopy 
(AFM), molecular dynamics (MD) simulations and 
magnetic tweezers. Our AFM studies yield insight in the 
overall compaction properties of the two proteins. The 
MD simulations using X-ray co-crystal structures yield 
an ensemble of equilibrated 'solution conformations' of 
the complexes. This gives information on the DNA- 
bending angle induced and the rigidity of this bend. 
With magnetic tweezers, we investigate the mechanical 
properties of single protein-DNA complexes. Elaborate 
analysis of measured force-extension curves enables us 
to quantify the number of proteins bound on single 
DNA molecules and to get a detailed understanding on 
the force-dependency of protein binding and the flexibility 
of the bending angle. 

MATERIALS AND METHODS 

Protein purification 

Crenl 

The Cren7 protein was purified from Escherichia coli 
strain BL21 (DE3) containing plasmid pET30a, including 
the gene encoding Cren7 (gene SSO6901) from 
S. solfataricus (19). Cells were grown in LB medium up 
to OD 600 x 0.4, and expression was induced using 0.5 mM 
IPTG at 37°C. Two hours after induction, cells were har- 
vested by centrifugation, washed with physiological saline 
[0.9% (w/v) NaCl] and resuspended in buffer A [50 mM 
Tris-HCl (pH 8.0), 10 mM ethylenediaminetetraacetic 
acid (EDTA), 10% glycerol, 10 mM P-mercaptoethanol]. 
Cells were lysed by sonication, and the cell lysate was 
centrifuged at 37 000 r. p.m. for 30 min at 4°C. The super- 
natant was applied to an SP column (GE Healthcare) and 
was eluted with buffer A containing 200 mM NaCl. The 
eluted protein was heated at 65°C for 30 min and 
centrifuged at 37 000 r. p.m. for 15 min. The supernatant 
was diluted to 100 mM NaCl in buffer A, loaded on a 
heparin column (GE Healthcare) and eluted with a 
linear gradient of 0.1-1 M NaCl in buffer A. The protein 
was eluted at -270 mM NaCl. 

Sul7 

A synthetic gene encoding the S. solfataricus Sul7 (Sso7d) 
protein (gene SSO10610) was constructed in expression 
vector pJexpress411 (supplied by DNA2.0). The DNA 
sequence used is available from the authors on request. 
The plasmid was transformed into E. coli strain BL21 
(DE3) Rosetta cells, grown in LB medium, and expression 
was induced using 0.2 mM IPTG at 37°C for 3 h. After cell 
lysis by sonication, the Sul7 protein was loaded onto a 
heparin column (GE Healthcare) and eluted with a 
linear gradient of 0-1 M NaCl in buffer B [20 mM 
Tris-HCl (pH 7.4), 1 mM dithiothreitol, 1 mM EDTA]. 
Fractions containing the Sul7 protein were pooled, 
concentrated and applied to a Superdex 200 gel filtration 
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column (GE Healthcare) and were eluted by isocratic flow 
with buffer B supplemented with 200 mM NaCl. 

Cren7 and Sul7 proteins were dialysed at 4°C against 
a storage buffer [20 mM HEPES (pH 7.5), 100 mM NaCl, 
10% glycerol, 10 mM P-mercaptoethanol], and they were 
stored at — 80°C until required. Protein concentrations 
were determined by ultraviolet absorbance at 280 nm, 
using a molar extinction coefficient £cv«7 = 8250 M -1 cm -1 
(19) and e Sul7 = 8300 M -1 cm -1 (21). 

DNA constructs 
AFM 

pUC18 was propagated in E. coli strain XL 10 and was 
purified (Qiagen plasmid midi kit). Digestion with 
EcoRI and XmnI resulted in a mixture of 845- and 
1841-bp linear fragments. 

Magnetic tweezers 

pBluescript II KS + (Stratagene) was digested with Sail 
and Hindlll enzymes to generate a 2947-bp linear 
fragment. The cohesive end at the Sail digestion site was 
labelled with dUTP-digoxigenin (DIG), using the Klenow 
fragment (Fermentas). Subsequently, the fragments were 
purified, labelled with dUTP-biotin at the extremity, 
resulting from Hindlll digestion using dATP, dCTP and 
dGTP and the Klenow fragment and were purified again. 

MD simulations 

MD simulations of 100 ns were performed in an explicit 
solvent environment as previously described (34) using 
initial coordinates from the co-crystal structures of 
Cren7-DNA and Sul7-DNA (PDB codes 3LWH and 
1BNZ, respectively), and an ideal B-DNA duplex with 
identical sequence was used for the bare DNA simulation. 
A detailed description can be found in the Supplementary 
'Materials and Methods' section (see Supplementary 
Data). 

AFM experiments 

Freshly cleaved mica was incubated with 10 (il of 0.05% 
(w/v) poly-L-lysine, rinsed with MilliQ water and dried 
with nitrogen gas. Protein-DNA complexes were formed 
by incubating 50 ng of DNA with varying concentrations 
of protein in 10 ul buffer containing lOmM HEPES 
(pH 7.5) and lOOmM NaCl for lOmin at room temperature 
(~23°C). After incubation, this mixture was diluted 5-fold 
in water and was directly deposited onto poly-L-lysine 
coated mica, rinsed with MilliQ water and dried with 
nitrogen gas. Images were collected on a NanoScopelll 
AFM (Digital Instruments, Santa Barbara, CA) using 
tapping mode in air (micro cantilevers, Olympus 
MCL-AC240TS-W2, resonance frequency 70 kHz, spring 
constant 2N/m) at a frequency of 2 Hz and were flattened 
using Nanoscope software (Veeco Instruments). 

Magnetic tweezers experiments 

Magnetic tweezers experiments were performed on a 
custom-built instrument described earlier (35). Images 
were acquired with a CCD camera (Pulnix TM-6710CL) 
at 60 Hz, and real time image processing was done using 



custom-developed LabView software (National 
Instruments). Magnets were positioned in the optical 
axis of the microscope and were controlled by a stepper 
motor-based translational stage (M-126, Physik 
Instrumente). Force was calibrated by measuring fluctu- 
ations of a tethered DNA-bead construct and was 
calculated according to the equipartition theorem (36). 

The bottom slide of the flow cell was pre-coated with 1 % 
(w/v) polystyrene in toluene. Subsequently, the flow cell 
was filled with 1 (ig/ml anti- DIG antibodies (Roche) and 
was incubated for 2 h at 4°C. Next, the flow cell was flushed 
with 2% (w/v) bovine serum albumin (BSA) and 1% (v/v) 
Tween-20 solution and was incubated overnight at 4°C. 
The flow cell was flushed with buffer I [10 mM HEPES 
(pH 7.6), 100 mM KAc, 0.2% (w/v) BSA, 0.1% (v/v) 
Tween-20], filled with 20ng/ml DNA (functionalized with 
biotin and DIG) in buffer I and incubated for lOmin. One 
microlitre of streptavidin coated superparamagnetic beads 
with a diameter of 1 |nn (DYNAL Myone) was washed 
three times in lxTE [10 mM Tris-HCl (pH 8.0), 1 mM 
EDTA], resuspended in 500 (il buffer I, flushed into the 
flow cell and incubated for 10 min. Before measurements, 
the flow cell was washed thoroughly (~10 times the volume 
of a flow cell) with buffer II [10 mM HEPES (pH 7.5), 
0.2% (w/v) BSA, 0.1% (v/v) Tween-20] with addition of 
either 25 or lOOmM NaCl (with or without lOmM MgCl 2 
added). Cren7 and Sul7 stock solutions were diluted in 
buffer II to the desired protein concentration 
(40-800 nM), and they were incubated for 10 min before 
starting measurements. Force-extension curves of single 
double-stranded DNA (dsDNA) molecules and protein- 
DNA complexes were measured by increasing the force to 
3.5 pN during 30 s and subsequently reducing the force to 
zero in 10 s. All measurements were performed at room 
temperature (~23°C). 

Data analysis 

Force-extension curves of single dsDNA molecules and 
protein-DNA complexes were quantified using the 
worm-like chain model (WLC). This model describes the 
mechanical properties of DNA and is defined by its intrin- 
sic properties: the persistence length L P and the contour 
length L 0 (37): 

\L P J\4(\ - z/L 0 y 4 LqJ 

Here, F denotes the applied force, z the end-to-end 
distance of the DNA molecule, k B the Boltzmann 
constant and T the absolute temperature. The effect of 
protein-DNA interactions on the DNA can be quantified 
by the fitting parameters L P and L 0 from F—z curves of 
protein-DNA complexes. 

The number of proteins that is bound to the DNA can 
be directly related to the change in end-to-end distance of 
the DNA molecule. Intercalation of a protein in the minor 
groove of the DNA that induces a bend in the DNA 
affects the end-to-end distance with a change in length 
Az (Figure 1). This change in end-to-end distance is 
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related 
X{F) = . 

Az = - 



to the force-dependent 
as follows (38,39): 



1 



A. • C(a) 



deflection length 



(2) 



where L P is the persistence length of bare DNA, and C(a) 
is a geometric coefficient depending on the bending angle 
a: 



C(a) = 8 



['-«©] 



(3) 



Implementing the reduction of length Az [Equation (2)] 
into the WLC [Equation (1)], and noting that 
1/(4(1 — z/Lq) 2 ) ^> z/Lo — 1/4 in our force regime, leads 
to a relation where the end-to-end distance z is directly 
related to the average number of proteins bound N B (38): 



z(F,N B ) 



1 



l_X{F) 
2 L P 



l+N B -C(a)- 



(4) 



This leads to an expression for the apparent persistence 
length, normalized by the number of proteins bound: 



l7(n b ) 



(\+N B .%C(a)) 



(5) 



The number of proteins bound along single DNA mol- 
ecules can now be determined by a model describing the 
mechanical response of DNA with multiple independent 
bends [Equation (4)]. It should be noted that this model 
only holds when the DNA is sufficiently stretched, such 
that the average distance between bound proteins is larger 
than the force-dependent deflection length, that is, 
N B < ^. In this force regime, bound proteins are distant 
enough to prevent the induced bends from interacting with 
each other. The aforementioned condition follows from an 
exact calculation that we performed for arbitrary linker 
lengths for a special case of protein arrangement (zigzag 
geometry). As for shorter linker lengths, the end-to-end 
distance depends on the rotational setting of adjacent 
proteins, our aforementioned theory can lead to under- 
or overestimation of the number of bound proteins 
outside the allowed regime. 



protein-DNA complexes using AFM. As can be seen in 
Figure 2, addition of either Cren7 or Sul7 results in DNA 
conformations that are more compact than bare DNA 
molecules. An increasing amount of protein results in a 
higher degree of compaction. It is likely that the observed 
compaction is caused by bends induced by the proteins. 
However, the small size of the proteins (~7 kDa) limits the 
unambiguous identification and quantitative analysis of 
individual proteins on the DNA. To gain more insights 
into the structure and dynamics of the bends induced by 
the two proteins, we carried out MD simulations of 
protein-DNA complexes in solution. 

Bending angle determination by MD simulations 

The co-crystal structures of Cren7-DNA and Sul7-DNA 
complexes showed that both proteins induce a bending in 
the DNA of ~50-60°. However, a complex in a crystal is 
fixed in a single conformation and may reflect a rare con- 
formation favoured by crystal packing. To evaluate 
whether the bending angles observed in the co-crystal struc- 
tures reflect the average conformation in solution, we per- 
formed MD simulations in solution using the co-crystal 
structures as initial coordinates. The MD simulations 
yielded stable structures except for a small deviation in 
the Sul7-DNA complex at the end of the simulation run 
because of rigid body motions of the C-terminal region 
(see Supplementary Figure SI for the root mean square 
deviations). The overall flexibility of the DNA in all three 
systems examined using root mean square fluctuations was 
found to be similar. Additionally, the extent of DNA 
bending was investigated by calculating the roll angle 
(40). Figure 3 shows the probability distributions of the 
roll angles during the simulation. Fitting the distributions 
to a Gaussian, we find that both protein-DNA complexes 
exhibit similar bending angles of a Cr en7 = 47.5 ± 5.6° 
and ocsun = 45.3 ± 6.7°, which slightly deviate from the 
bending angle found in the co-crystal structures. It is 
important to note that the width of the roll angle distribu- 
tions (which reflects the flexibility of the complex) of both 
protein-DNA complexes is small and comparable with the 
width of the roll angle distribution of bare DNA 

(^bareDNA = 7.3 ± 5.7°). 



RESULTS AND DISCUSSION 

To investigate the architectural properties of Cren7 and 
Sul7, we used a combination of different single-molecule 
techniques. We used AFM to visualize individual protein- 
DNA complexes, MD simulations to evaluate the bending 
angle of the protein-DNA complexes in solution and 
magnetic tweezers experiments to quantify the physical 
properties of single protein-DNA complexes and their 
force dependency. 

Effect of Cren7 and Sul7 on the conformation of 
single DNA molecules 

To investigate the effect of Cren7 and Sul7 on the con- 
formation of single DNA molecules, we visualized 
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Figure 2. AFM images of bare DNA, Cren7-DNA complexes and Sul7- 
DNA complexes. The protein-DNA complexes are more compact than 
the bare DNA molecules. Increasing the protein concentration results in 
more DNA compaction. The DNA is an equimolar mixture of both 
1841- and 845-bp linear fragments. Images are 500 x 500 nm in size. 
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Figure 3. Probability distribution of the roll angle of the kinked base 
pair step in the bare DNA ('blue'), the DNA in the Cren7-DNA 
complex ('green') and the DNA in the Sul7-DNA complex ('red'). 
Distributions are obtained from the final 90 ns of the MD simulations 
and are fitted to a Gaussian (ai, areDNA = 7.3 ± 5.7°, a Cren7 = 
47.5 ± 5.6° and a Su n = 45.3 ± 6.7°). Insets: the average conformations 
of the bare DNA, the Cren7-DNA complexes and Sul7-DNA 
complexes are shown and aligned with two extreme conformations to 
illustrate the flexibility of the complexes. Images are generated using 
VMD software (22). 



Quantification of DNA compaction by Cren7 and Sul7 

To quantify the molecular mechanisms underlying the 
observed DNA compaction, we carried out single- 
molecule micromanipulation experiments. In these experi- 
ments, we measured the force-extension relation of single 
dsDNA molecules over a range of concentrations for both 
Cren7 and Sul7 (0, 40, 80, 200, 400 and 800 nM) (Figure 
4A and B). In both cases, the addition of protein results in 
a reduction of z (compaction), which decreases progres- 
sively on raising the protein concentration. No significant 
difference is observed between Cren7 and Sul7 as they 
compact the DNA to a similar extent, which confirms 
the compaction observed in our AFM studies. The reduc- 
tion of the end-to-end distance (z) saturates at a protein 
concentration of ~400nM (Figure 4C and D). At the 
highest concentration (800 nM), the end-to-end distance 
increases slightly, which could point to interactions 
between adjacent proteins, causing a mild extension of 
the compacted DNA molecules. However, incubation 
with higher protein concentrations (up to 1600 nM) did 
not induce a further increase in end-to-end distance 
('data not shown'). 

To quantify the effects of the protein on the mechanical 
properties of DNA, the WLC model [Equation (1)] was 
fitted to the force-extension curves. Although the contour 
length stays constant (L 0 = 1.0 ± 0.05 Lim), the apparent 



persistence length (38,39) decreases drastically 
(Figure 5A). At high concentrations, the apparent persist- 
ence length reaches a minimum of L P = 8.8 ± 1.3 nm for 
Cren7 and L P = 9.0 ± 1.6 nm for Sul7, which is more than 
four times smaller than the measured persistence length of 
bare DNA, L P = 40.5 ± 11.2nm. Thus, binding of these 
proteins results in compaction of DNA by reducing the 
apparent persistence length. 

It is interesting to note that higher protein concentra- 
tions do not result in stiffening of the DNA. Previous 
studies of similar DNA-bending proteins (e.g. the bacter- 
ial proteins HU and IHF and eukaryotic HMG proteins) 
have shown that these proteins compact DNA by forming 
bends at low protein concentrations (41-44). However, 
HU and HMG proteins enter a second binding regime 
at high concentrations, forming stiff filaments because of 
close adjacent binding. It has been suggested that IHF 
behaves similarly (2) when binding non-specifically along 
DNA (44,45), but this has not yet been experimentally 
confirmed. The bimodal binding of HU, causing a transi- 
tion from a softening to a stiffening mode, has been shown 
to be salt dependent and was only observed at NaCl con- 
centrations <150mM (46). As such salt-dependent effects 
might exist also for Cren7 and Sul7, we measured 
force-extension curves under different ionic conditions. 
Previous studies have already shown that the binding 
affinity of Sul7 decreases with increasing ionic strength 
(47). Measurements in a buffer containing a lower salt 
concentration (25 mM NaCl) showed that the decrease 
in apparent persistence length occurs at much lower 
protein concentration (Figure 5B), which suggests an 
increased binding affinity, because of enhanced electro- 
static protein-DNA interactions at low ionic strength, as 
expected. Again at the highest protein concentration, the 
persistence length did not increase, and a stiffening mode 
was not observed. These findings suggest that Cren7 and 
Sul7 do not stiffen DNA in a second binding mode, even 
at protein concentrations where high DNA coverage is 
expected. It has been shown for bacterial H-NS-like 
proteins that these proteins exhibit two binding modes 
(48-50), which can be switched by the addition of 
MgCl 2 (51,52). To test whether divalent ions have an 
effect on the mode of binding of Cren7 and Sul7, we per- 
formed measurements in the presence of MgCl 2 . 
These experiments did not reveal an effect on the 
binding behaviour of Cren7 and Sul7, other than a 
decrease in binding affinity (see Supplementary Figure 
S2). This reinforces that the binding is dominated by elec- 
trostatic interactions, which are reduced by increasing the 
ionic strength. 

Quantification of the number of proteins bound 

The number of proteins bound to the DNA within the 
range of measured force and protein concentrations is 
determined by applying Equation (4) to the measured 
force-extension curves (Figure 4A and B) using the 
bending angle a determined from the MD simulations 
(a c , M7 = 47.5° and a Sul7 = 45.3, see Figure 3). As 
expected, the number of proteins bound on the DNA 
(N B ) increases with increasing protein concentration 
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Figure 4. Force-extension curves of protein-DNA complexes. Each data point represent the average of a number of DNA molecules measured in 
buffer containing 100 mM NaCl (n = 26 for bare DNA molecules and n = 6-14 for DNA-protein complexes). Data are fitted with the WLC model 
[Equation (1)]. Error bars represent the standard deviation. (A and B) Averaged force-extension curves of DNA at different protein concentrations of 
Cren7 and Sul7. (C and D) End-to-end distance as a function of bulk protein concentration represented at constant forces. Progressive addition of 
protein results in a decrease of the end-to-end distance z, and thus compaction of DNA molecules. 



(see Figure 6). Considering the condition for which our 
model is valid, Nb < ^ (see 'Materials and Methods' 
section), we note that this condition holds for the 40 and 
80 nM measurements, where proteins are still bound 
distant enough from each other. For the lowest concen- 
trations, the number of bound proteins is found to be 
constant throughout the measured force range (Figure 
6), from which we conclude that forces up to 3.5 pN do 
not affect the number of proteins on the DNA. At higher 
protein concentrations, the high occupancy may introduce 
a systematic error in the calculated N B , but as in the case 
of lower concentrations, we still find a force independent 
value for N B . This is a surprising observation, as one 
would expect a significant decrease in binding affinity 
with applied force. The energy it costs to form a bend is 
force-dependent, and for the highest force (3.5 pN), this 
would be E k!nk = ~k B T {see Supplementary Data). As the 



dissociation constant K,/is related to E kink , and the binding 
free energy AG 0 as follows: 



~ eXP ( k B T 



(6) 



One would expect the dissociation constant to increase 
with a factor 2.7 for both Cren7 and Sul7 when the force 
is increased from 0 to 3.5 pN, resulting in a lower protein 
occupancy at higher forces. Apparently, these proteins 
bind stably enough to withstand such forces on 
our experimental timescale (~30s). Note that measure- 
ments were reversible for all measured protein concentra- 
tions, which could be explained by off-rates much larger 
than the experimental timeframe (see Supplementary 
Figure S3). Even holding the protein-DNA complexes 
at 3.5 pN for 10 min did not result in an increase in 
z or a different retraction curve, showing that the 
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Figure 6. Number of bound proteins (N B ) on a single DNA molecule (2947 bp) as a function of force, calculated from the force-extension curves 
according to Equation (4). The number of proteins increases with increasing protein concentration, but is unaffected by forces up to 3.5 pN. Lines 
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proteins stay stably bound even on longer timescales 
(see Supplementary Figure S4). Unexpectedly, low off- 
rates of proteins bound to DNA have also been observed 
in previous studies on HU and HMG protein-DNA 
complexes in protein-free solution (46,53). Apparently, 
there is a large energy barrier for the proteins to dissociate, 
which is not reached by applying a force of 3.5 pN. Such 
strong binding could be crucial for maintaining the integ- 
rity of the DNA in extreme environmental conditions of 
the organism's natural habitat. Cellular processes that 
apply local tension on the DNA, such as transcription, 
replication and repair, act on timescales that are much 
shorter compared with our experimental timescale. This 
suggests that these proteins are relatively stably bound 
to the DNA in vivo. However, cellular machineries 
involved in replication and transcription, such as helicases 
and polymerases, may directly drive the proteins to dissoci- 
ate (54), allowing access to the DNA track. 

Determining the number of proteins bound to DNA 
and how this is affected by force on the DNA has been 



a topic of interest in various recent studies (55-58). It is 
important in the light of understanding the contribution 
of the individual bound proteins to overall genome organ- 
ization and how forces induced by cellular processes affect 
this. Theoretical models have described the number of 
proteins bound to single DNA molecules in relation to 
applied force. These models showed that the force 
response of the bound proteins is highly dependent on 
the flexibility of the protein-DNA complex if the 
binding is relatively strong (55). With the use of thermo- 
dynamic Maxwell relations (57), the change in the number 
of bound proteins can be determined, given the condition 
that the protein-DNA interactions are in thermodynamic 
equilibrium. Analysis with the thermodynamic Maxwell 
relations of single-DNA force-extension data with 
DNA-bending proteins from E. coli, HU and FIS, 
suggested that these proteins are driven off the DNA by 
an applied force, even at relatively small forces <1 pN 
(56). It should be noted, however, that this approach 
does not take into account that the degree of protein 
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induced DNA bending could be altered by DNA tension. 
Applying force to the protein-DNA complexes could 
reduce the degree of bending, which would result in a 
net mechanical response, similar to that of protein unbind- 
ing (41). Thus, the proposed unbinding of proteins in 
those studies might actually be (in part) caused by a 
decrease of bending angle rather than protein unbinding. 

Do Cren7 and Sul7 induce flexible or static bends? 

Previous studies of bacterial chromatin protein HU have 
shown that bends of HU-DNA complexes are highly 
flexible (59) [varying from 0 to 180° (41)]. Although a 
similar flexible hinge model was initially proposed for eu- 
karyotic HMG domains (42), further experiments showed 
a more narrow angle distribution of high mobility group B 
(HMGB) protein-DNA complexes (31). Non-specific 
binding of IHF may induce flexible bends analogous to 
HU, whereas the site-specific binding of IHF is stabilized 
by interaction between the DNA flanking the kink sites 
and the body of the protein, making this bend static (32). 
Based on the analogy of Cren7 and Sul7 with those 
non-specific DNA-bending proteins, one would anticipate 
that these proteins induce bends with a certain degree of 
flexibility. Applying a force on such flexible bends would 
decrease the bending angle. Equation (4) yields an increase 
of the extension of the DNA for both a reduction in the 
bending angle and a reduction of N B , which are indistin- 
guishable from each other. Surprisingly, the number of 
bound proteins N B is unaffected by the applied force, 
which implies that the bending angles are independent of 
forces <3.5 pN. The MD simulations of the protein-DNA 
complexes in solution at zero force show that the distri- 
butions of roll angles are narrow (Figure 3) 
(oLcrenJ = 47.5 ± 5.6° and a Su/7 = 45.3 ± 6.7°) and com- 
parable in width with the roll angle distribution of bare 
DNA (ahareDNA = 7.3 ± 5.7°). Both observations show 
that Cren7 and Sul7 bind with a rigid bend, in contrast 
with flexible hinges induced by analogous DNA-bending 
proteins. If a flexible angle is essential in permitting tight 
packing of proteins, this could explain the absence of a 
second-binding mode in which the DNA is stiffened. 

CONCLUSIONS 

The protein-DNA interactions of Cren7 and Sul7 are 
similar. Like in previous biochemical assays, we found 
that they bind with comparable affinity. Here, we have 
demonstrated that they induce a similar degree of com- 
paction of DNA, bind independent of forces up to 3.5 pN 
and bend DNA with a non-flexible bend. Whether the 
bending and binding of Cren7 and Sul7 also resists 
higher forces remains to be seen and could be assessed 
with other single-molecule techniques that permit 
DNA-stretching at higher forces. Visualization of fluores- 
cent proteins in combination with such techniques (60) 
could confirm our findings in the low force regime. 
Although no difference has yet been observed, given the 
fact that Cren7 and Sul7 coexist in the same organisms, 
one would still expect them to have a functional difference 
in vivo. As mentioned earlier, Cren7 and Sul7 can both be 



methylated at several lysine residues. Such post- 
translational modifications of nucleoid-associated 
proteins can play an important role in regulation of inter- 
actions between individual proteins on the DNA, as well 
as the interaction with other proteins, such as transcrip- 
tion factors or other chromatin proteins. How this affects 
the role of Cren7 and Sul7 in chromatin organization, and 
possibly gene regulation, remains to be investigated. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Figures 1^1, Supplementary Information, 
Supplementary Materials and Methods and 
Supplementary References [61-71]. 

ACKNOWLEDGEMENTS 

The authors thank Dr Nora Goosen and Geri Moolenaar 
for technical assistance and helpful discussions. They 
thank Federica Galli and Tjerk Oosterkamp for their 
maintenance and supervision of the Bio-AFM Lab. 
Biljana Petrovic-Stojanovska is acknowledged for tech- 
nical assistance with the purification of Sul7. They thank 
Dr Li Huang (Chinese Academy of Sciences, Beijing, 
China) for kindly providing the expression vector 
pET30a, including the gene encoding Cren7. 

FUNDING 

Netherlands Organization for Scientific Research (NWO) 
[864.08.001 to R.T.D.]; Federation of European 
Biochemical Societies (to R.S.); Human Frontiers 
Science Program (to J.v.N.); Council of Scientific and 
Industrial Research (CSIR), India (to G.S.). Funding for 
open access charge: NWO. 

Conflict of interest statement. None declared. 



REFERENCES 

1. Luijsterburg,M.S., White,M.F., van Driel.R. and Dame.R.T. 
(2008) The major architects of chromatin: architectural proteins 
in bacteria, archaea and eukaryotes. Crit. Rev. Biochem. Mot. 
Biol, 43, 393^18. 

2. Dame,R.T. (2005) The role of nucleoid-associated proteins in the 
organization and compaction of bacterial chromatin. Mol. 
Microbiol, 56, 858-870. 

3. Luijsterburg,M.S., Noom,M.C, Wuite,G.J. and Dame,R.T. (2006) 
The architectural role of nucleoid-associated proteins in the 
organization of bacterial chromatin: a molecular perspective. 

J. Struct. Biol, 156, 262-272. 

4. Bell,S.D. and White,M.F. (2010) Archaeal Chromatin 
Organization. In: Dame.R.T. and Dorman,C.J. (eds), Bacterial 
Chromatin. Springer, Dordrecht, the Netherlands, pp. 205-217. 

5. Sandman, K. and ReeveJ.N. (2005) Archaeal chromatin proteins: 
different structures but common function? Curr. Opin. Microbiol, 
8, 656-661. 

6. ReeveJ.N., Bailey, K.A., Li,W.-t., Marc,F., Sandman,K. and 
Soares,D.J. (2004) Archaeal histones: structures, stability 
and DNA binding. Biochem. Soc. Trans., 32, 227-230. 

7. Pereira,S.L., Grayling,R.A., Lurz.R. and ReeveJ.N. (1997) 
Archaeal nucleosomes. Proc. Natl Acad. Sci. USA, 94, 
12633-12637. 



204 Nucleic Acids Research, 2013, Vol. 41, No. 1 



8. Cubonova,L., Sandman,K., Hallam,S.J., Delong,E.F. and 
ReeveJ.N. (2005) Histones in crenarchaea. /. BacterioL, 187, 
5482-5485. 

9. Driessen.R.P.C. and Dame.R.T. (2011) Nucleoid-associated 
proteins in Crenarchaea. Biochem. Soc. Trans., 39, 116-121. 

10. Wardleworth,B.N., Russell,R.J., Bell,S.D., Taylor.G.L. and 
White,M.F. (2002) Structure of Alba: an archaeal 
chromatin protein modulated by acetylation. EMBO J., 21, 
4654-4662. 

11. Lurz,R., Grote,M., DijkJ., Reinhardt,R. and Dobrinski,B. (1986) 
Electron microscopic study of DNA complexes with proteins 
from the Archaebacterium Sulfolobus acidocaldarius. EMBO J., 5, 
3715-3721. 

12. Edmondson.S.P., Kahsai,M.A., Gupta.R. and ShriverJ.W. (2004) 
Characterization of Sac 10a, a hyperthermophile DNA-binding 
protein from Sulfolobus acidocaldarius. Biochemistry, 43, 
13026-13036. 

13. Kahsai,M.A., Vogler,B., Clark,A.T., Edmondson.S.P. and 
ShriverJ.W. (2005) Solution structure, stability, and flexibility 
of SsolOa: a hyperthermophile coiled-coil DNA-binding protein. 
Biochemistry, 44, 2822-2832. 

14. Napoli,A., Zivanovic,Y., Bocs,C, Buhler,C, Rossi,M., Forterre,P. 
and Ciaramella,M. (2002) DNA bending, compaction and 
negative supercoiling by the architectural protein Sso7d of 
Sulfolobus solfataricus. Nucleic Acids Res., 30, 2656-2662. 

15. Edmondson,S.P. and ShriverJ.W. (2001) DNA binding proteins 
Sac7d and Sso7d from Sulfolobus. Methods Enzymol., 334, 
129-145. 

16. McAfeeJ.G., Edmondson,S.P., Datta,P.K., ShriverJ.W. and 
Gupta,R. (1995) Gene cloning, expression, and characterization 
of the Sac7 proteins from the hyperthermophile Sulfolobus 
acidocaldarius. Biochemistry, 34, 10063-10077. 

17. McAfeeJ.G., Edmondson,S.P., Zegar,I. and ShriverJ.W. (1996) 
Equilibrium DNA binding of Sac7d protein from the 
hyperthermophile Sulfolobus acidocaldarius: fluorescence and 
circular dichroism studies. Biochemistry, 35, 4034-4045. 

18. Choli,T., Henning,P., Wittmann-Liebold.B. and Reinhardt,R. 
(1988) Isolation, characterization and microsequence analysis of 
a small basic methylated DNA-binding protein from the 
Archaebacterium, Sulfolobus solfataricus. Biochim. Biophys. Acta, 
950, 193-203. 

19. Guo,L., Feng,Y., Zhang.Z., Yao,H., Luo,Y., WangJ. and 
Huang,L. (2008) Biochemical and structural characterization of 
Cren7, a novel chromatin protein conserved among Crenarchaea. 
Nucleic Acids Res., 36, 1129-1137. 

20. Mai.V.Q., Chen,X., Hong,R. and Huang,L. (1998) Small 
abundant DNA binding proteins from the thermoacidophilic 
archaeon Sulfolobus shibatae constrain negative DNA supercoils. 
J. BacterioL, 180, 2560-2563. 

21. Baumann,H., Knapp.S., Lundback,T., Ladenstein.R. and Hard,T. 
(1994) Solution structure and DNA-binding properties of a 
thermostable protein from the archaeon Sulfolobus solfataricus. 
Nat. Struct. Biol, 1, 808-819. 

22. Humphrey,W., Dalke,A. and Schulten,K. (1996) VMD: visual 
molecular dynamics. J. Mol. Graph., 14, 33-38. 

23. Feng,Y., Yao,H. and WangJ. (2010) Crystal structure of the 
crenarchaeal conserved chromatin protein Cren7 and 
double-stranded DNA complex. Protein Set, 19, 1253-1257. 

24. Zhang.Z., Gong,Y., Guo,L., Jiang,T. and Huang,L. (2010) 
Structural insights into the interaction of the crenarchaeal 
chromatin protein Cren7 with DNA. Mol. Microbiol., 76, 
749-759. 

25. Gao,Y.G., Su,S.Y., Robinson.H., Padmanabhan,S., Lim,L., 
McCrary,B.S., Edmondson.S.P., ShriverJ.W. and Wang,A.H. 
(1998) The crystal structure of the hyperthermophile 
chromosomal protein Sso7d bound to DNA. Nat. Struct. Biol., 5, 
782-786. 

26. Robinson,H., Gao,Y.G., McCrary,B.S., Edmondson,S.P., 
ShriverJ.W. and Wang,A.H. (1998) The hyperthermophile 
chromosomal protein Sac7d sharply kinks DNA. Nature, 392, 
202-205. 

27. Agback,P., Baumann,H., Knapp,S., Ladenstein,R. and Hard,T. 
(1998) Architecture of nonspecific protein-DNA interactions in 
the Sso7d-DNA complex. Nat. Struct. Bio/., 5, 579-584. 



28. WurtzekO., Sapra,R., Chen,F., Zhu,Y., Simmons.B.A. and 
Sorek,R. (2010) A single-base resolution map of an archaeal 
transcriptome. Genome Res., 20, 133-141. 

29. Bell.S.D., Botting,C.H., Wardleworth,B.N., Jackson,S.P. and 
White,M.F. (2002) The interaction of Alba, a conserved archaeal 
chromatin protein, with Sir2 and its regulation by acetylation. 
Science, 296, 148-151. 

30. ThomasJ.O. and Travers,A.A. (2001) HMG1 and 2, and related 
'architectural' DNA-binding proteins. Trends Biochem. Sci., 26, 
167-174. 

31. Zhang,J., McCauley,M.J., Maher,L.J. 3rd, Williams, M.C. and 
Israeloff,N.E. (2009) Mechanism of DNA flexibility enhancement 
by HMGB proteins. Nucleic Acids Res., 37, 1107-1114. 

32. Swinger.K.K. and Rice,P.A. (2004) IHF and HU: 

flexible architects of bent DNA. Curr. Opin. Struct. Biol., 14, 
28-35. 

33. Bessa,E., Wintraecken,K., Geerling,A. and de Vries,R. (2007) 
Synergy of DNA-bending nucleoid proteins and macromolecular 
crowding in condensing DNA. Biophys. Rev. Lett., 2, 1-7. 

34. Priyakumar,U.D., Harika,G. and Suresh,G. (2010) Molecular 
simulations on the thermal stabilization of DNA by 
hyperthermophilic chromatin protein Sac7d, and associated 
conformational transitions. /. Phys. Chem. B, 114, 16548-16557. 

35. Kruithof,M., Chien,F., de Jager,M. and van NoortJ. (2008) 
Subpiconewton dynamic force spectroscopy using magnetic 
tweezers. Biophys. J., 94, 2343-2348. 

36. Strick,T.R., AllemandJ.F., Bensimon,D., Bensimon,A. and 
Croquette,V. (1996) The elasticity of a single supercoiled DNA 
molecule. Science, 271, 1835-1837. 

37. MarkoJ.F. and Siggia,E.D. (1995) Stretching DNA. 
Macromolecules, 28, 8759-8770. 

38. Kulic.I.M., Mohrbach,H., Thaokar,R. and Schiessel,H. (2007) 
Equation of state of looped DNA. Phys. Rev. E Stat. Nonlin. 
Soft Matter Phys., 75, 011913. 

39. KulicJ.M., Mohrbach,H., Lobaskin,V., Thaokar.R. and 
Schiessel,H. (2005) Apparent persistence length renormalization 
of bent DNA. Phys. Rev. E Stat. Nonlin. Soft Matter Phys., 72, 
041905. 

40. Lavery,R., Moakher,M., MaddocksJ.H., Petkeviciute,D. and 
Zakrzewska,K. (2009) Conformational analysis of nucleic acids 
revisited: Curves+. Nucleic Acids Res., 37, 5917-5929. 

41. van NoortJ., Verbrugge,S., Goosen,N., Dekker.C. and 
Dame,R.T. (2004) Dual architectural roles of HU: formation 

of flexible hinges and rigid filaments. Proc. Natl Acad. Sci. USA, 
101, 6969-6974. 

42. McCauley,M., Hardwidge,P.R., Maher,L.J. 3rd and 
Williams,M.C. (2005) Dual binding modes for an HMG domain 
from human HMGB2 on DNA. Biophys. J., 89, 353-364. 

43. Skoko,D., Wong.B., Johnson,R.C. and MarkoJ.F. (2004) 
Micromechanical analysis of the binding of DNA-bending 
proteins HMGB1, NHP6A, and HU reveals their ability to form 
highly stable DNA-protein complexes. Biochemistry, 43, 
13867-13874. 

44. Ali.B.M., Amit,R., Braslavsky,I., Oppenheim.A.B., Gileadi,0. and 
StavansJ. (2001) Compaction of single DNA molecules induced 
by binding of integration host factor (IHF). Proc. Natl Acad. Sci. 
USA, 98, 10658-10663. 

45. HolbrookJA., Tsodikov,O.V., Saecker,R.M. and Record.M.T. Jr 
(2001) Specific and non-specific interactions of integration host 
factor with DNA: thermodynamic evidence for disruption of 
multiple IHF surface salt-bridges coupled to DNA binding. 

/. Mol. Biol., 310, 379^101. 

46. Xiao,B., Johnson,R.C. and MarkoJ.F. (2010) Modulation of 
HU-DNA interactions by salt concentration and applied force. 
Nucleic Acids Res., 38, 6176-6185. 

47. Lundback,T. and Hard,T. (1996) Salt dependence of the free 
energy, enthalpy, and entropy of nonsequence specific DNA 
binding. /. Phys. Chem., 100, 17690-17695. 

48. Dame,R.T., Wyman,C. and Goosen,N. (2000) H-NS mediated 
compaction of DNA visualised by atomic force microscopy. 
Nucleic Acids Res., 28, 3504-3510. 

49. Dame,R.T., Noom,M.C. and Wuite,G.J. (2006) Bacterial 
chromatin organization by H-NS protein unravelled using dual 
DNA manipulation. Nature, 444, 387-390. 



Nucleic Acids Research, 2013, Vol. 41, No. 1 205 



50. Wiggins,P.A., Dame,R.T., Noom.M.C. and Wuite.G.J. (2009) 
Protein-mediated molecular bridging: a key mechanism in 
biopolymer organization. Biophys, J., 97, 1997-2003. 

51. Liu,Y., Chen,H., Kenney,L-J- and YanJ. (2010) A divalent switch 
drives H-NS/DNA-binding conformations between stiffening and 
bridging modes. Genes Dev., 24, 339-344. 

52. Lim,C.J., Whang,Y.R., Kenney,L.J. and YanJ. (2012) Gene 
silencing H-NS paralogue StpA forms a rigid protein filament 
along DNA that blocks DNA accessibility. Nucleic Acids Res., 
40, 3316-3328. 

53. GrahamJ.S., Johnson.R.C. and MarkoJ.F. (2010) 
Concentration-dependent exchange accelerates turnover of 
proteins bound to double-stranded DNA. Nucleic Acids Res., 39, 
2249-2259. 

54. Bustamante,C, Cheng,W. and Mejia.Y.X. (2011) Revisiting the 
central dogma one molecule at a time. Cell, 144, 480—497. 

55. YanJ. and MarkoJ.F. (2003) Effects of DNA-distorting proteins 
on DNA elastic response. Phys. Rev. E Stat. Nonlin. Soft Matter 
Phys., 68, 011905. 

56. Xiao,B., Zhang,H., Johnson,R.C. and MarkoJ.F. (2011) 
Force-driven unbinding of proteins HU and Fis from DNA 
quantified using a thermodynamic Maxwell relation. Nucleic Acids 
Res., 39, 5568-5577. 

57. Zhang,H. and MarkoJ.F. (2008) Maxwell relations for 
single-DNA experiments: monitoring protein binding and 
double-helix torque with force-extension measurements. Phys. 
Rev. E Stat. Nonlin. Soft Matter Phys., 77, 031916. 

58. Liebesny.P., Goyal.S., Dunlap,D., Family, F. and Finzi,L. (2010) 
Determination of the number of proteins bound non-specifically 
to DNA. /. Phys. Condens. Matter, 22, 414104. 

59. Swinger,K.K., Lemberg,K.M., Zhang.Y. and Rice,P.A. (2003) 
Flexible DNA bending in HU-DNA cocrystal structures. EMBO 
J., 22, 3749-3760. 

60. Candelli,A., Wuite,G.J. and Peterman.E.J. (2011) Combining 
optical trapping, fluorescence microscopy and micro-fluidics for 
single molecule studies of DNA-protein interactions. Phys. Chem. 
Chem, Phys., 13, 7263-7272. 

61. Brooks,B.R., Brooks.C.L. 3rd, Mackerell,A.D. Jr, Nilsson,L., 
Petrella,R.J., Roux,B., Won,Y., Archontis,G, Bartels,C, 



Boresch.S. et al. (2009) CHARMM: the biomolecular simulation 
program. J. Comput. Chem., 30, 1545-1614. 

62. Foloppe.N. and MacKerellJ.A.D. (2000) All-atom empirical force 
field for nucleic acids: I. Parameter optimization based on small 
molecule and condensed phase macromolecular target data. 

/. Comput. Chem., 21, 86-104. 

63. MacKerell,A.D. and Banavali,N.K. (2000) All-atom empirical 
force field for nucleic acids: II. Application to molecular 
dynamics simulations of DNA and RNA in solution. /. Comput. 
Chem., 21, 105-120. 

64. MacKerell,A.D., Bashford.D., Dunbrack,R.L., EvanseckJ.D., 
Field.M.J., Fischer,S., GaoJ., Guo,H., Ha,S., Joseph- 
McCarthy^, et al. (1998) All-atom empirical potential for 
molecular modeling and dynamics studies of proteins. J. Phys. 
Chem. B, 102, 3586-3616. 

65. Mackerell.A.D. Jr, Feig,M. and Brooks,C.L. 3rd (2004) Extending 
the treatment of backbone energetics in protein force fields: 
limitations of gas-phase quantum mechanics in reproducing 
protein conformational distributions in molecular dynamics 
simulations. J. Comput. Chem., 25, 1400-1415. 

66. Jorgensen,W.L., ChandrasekharJ., MaduraJ.D., Impey,R.W. and 
Klein.M.L. (1983) Comparison of simple potential functions for 
simulating liquid water. J. Chem. Phys., 79, 926-935. 

67. RyckaertJ.P., Ciccotti,G. and Berendsen,H.J.C. (1977) Numerical 
integration of the Cartesian equations of motion of a system with 
constraints: molecular dynamics of n-alkanes. /. Comput. Phys., 
23, 327-341. 

68. Essmann,U., Perera,L., Berkowitz.M., Darden,T., Lee,H. and 
Pedersen,L. (1995) A smooth particle mesh Ewald method. 

/. Chem. Phys., 103, 8577-8593. 

69. Darden.T., Perera,L., Li,L. and Pedersen,L. (1999) New tricks for 
modelers from the crystallography toolkit: the particle mesh 
Ewald algorithm and its use in nucleic acid simulations. Structure, 
7, R55-R60. 

70. Hoover.W.G. (1985) Canonical dynamics: equilibrium phase-space 
distributions. Phys. Rev. A, 31, 1695-1697. 

71. Feller,S.E., Zhang,Y., Pastor,R.W. and Brooks,B.R. (1995) 
Constant pressure molecular dynamics simulation: the Langevin 
piston method. /. Chem. Phys., 103, 4613^1621. 



